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Overview 

>  What  is  reliability  ?? 

>  Basics  of  reliability  methods  for  repairable  and  non¬ 
repayable  systems 

>  Estimation  of  PDF  of  Time  Between  Failures  (TBF) 
using  limited,  censored  data 

•  “Frequentist”  approach  (Method  1) 

•  Bayesian  updating  approach  (Method  2) 

S  “Enhances”  data  with  expert  opinion 
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What  is  Reliability? 


Reliability  at  time  t  is  the  probability  that  the  system 

has  not  failed  before  time  t. 
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failure 


— ► 

time 


R(t)  =  P(T>t)  =  l-P(T<t) 
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Reliability  Basics  for 
Non-Repairable 

Systems 
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Reliability  of  Non-Repairable  Systems 


failure 


o 


t 


time 


R(t)  =  P(T>t)  =  l-P(T<t)=>R(t)  =  l-F(t) 


(1) 


/ 

Failure 

Rate 


P(t  <T  <t  +  dt/T  >t)  P{t<T<t  +  dt ) 


dt 


dt*p(r>t) 


_  F(t  +  dt)-F(t) 
dt*R(t) 


m =4j 

R(t) 


(2) 
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Reliability  of  Non-Repairable  Systems 


failure 
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time 


dR 

dt 


due  to  (2) 
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=  —Adt  =>  c/(ln  R)  =  —Adt  =>  In 


rm' 
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—  f  Adt 


l 

R(t)=  exp[-f  Adt] 
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All  we  need  is  the 
failure  rate 
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Reliability  of  Non-Repairable  Systems 
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Reliability  of  Non-Repairable  Systems 
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Reliability  of  Non-Repairable  Systems 

j= 1 

i- 1 

j= i 

N<, 

R.  —  e  ' 
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Reliability  Calculation 

All  we  need  for  calculating  the 
reliability  of  a  system  (non-repairable 
or  repairable)  is  the  system  PDF  of 
time  to  failure  (TTF) 

We  use  : 

>  Data  to  estimate  the  PDF  of  TTF  for  each  component 

>  Monte  Carlo  simulation  to  estimate  the  PDF  of  TTF 
for  the  system 
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Basics  of  Reliability  Methods 
(Repairable  Systems) 
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Reliability  of  Repairable  Systems 
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Reliability  of  Repairable  Systems 


kA t  system  failures  in  At 
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At  segment 
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At 
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:  Probability  of  failure  in  At  (very  small  if  At  — »  0  ) 


/ 

#  of  systems  (vehicles)  in  fleet 


If  p  is  independent  of  At  segment 


If p  depends  on  At  segment  t 


> 


Homogeneous  Poisson 
Process  (HPP) 

Non-Homogeneous 
Poisson  Process  (NHPP) 
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Homogeneous  Poisson  Process  (HPP) 


pn  :  Expected  (average)  #  of  failures  /  system  in  (0,  T] 

Probability  of 
failure  in  At 


Define  : 


AT  =  pn 


A  :  Average  #  of  failures  per  system  per  unit  time 

( failure  rate,  hazard  rate,  intensity  rate, 
repair  rate ) 
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Homogeneous  Poisson  Process  (HPP) 

Reliability  Calculation  for  HPP  at  T  =  11  At 


♦>  Assume  statistically  independent  events  at  each  At 


r(t)  =  (1- p\l- p)=(l- p)" 
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Non-Homogeneous  Poisson  Process  (NHPP) 

Reliability  Calculation  for  NHPP  at  T  =  ft  At 
❖  Assume  statistically  independent  events  at  each  At 


=  (l  —  A^At  )(l  —  X-,  At )  •  •  •  (l  —  AnAt) 


£) 


Same  formula  with  non¬ 
rep  air  able  systems 


if  At  — ^  0 


Cumulative  Hazard 
Rate 
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Reliability  Calculation 

Series  Systems 


R  =  R,  * R-,  *■■■* R  = 
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Calculation  of  Hazard  Rate  A (t) 

Mf)\ 


A 


n  T  time 

At 

Determine  average  #  of  failures  kAt  (t )  among  N  systems 
in  each  At .  Then  : 
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Calculation  of  Hazard  Rate  X(t) 


For  constant  hazard  rate  systems  : 


MTBF 


For  most  engineering  systems,  the  hazard  rate  IS  NOT 
constant.  To  estimate  it,  we  need  the  PDF  of  TTF  for 

each  component. 
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Reliability  Calculation 

All  we  need  for  calculating  the 
reliability  of  a  system  (non-repairable 
or  repairable)  is  the  system  PDF  of 
time  to  failure  (TTF) 

We  use  : 

>  Data  to  estimate  the  PDF  of  TTF  for  each  component 

>  Monte  Carlo  simulation  to  estimate  the  PDF  of  TTF 
for  the  system 
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Estimation  of  the  PDF  (or 
CDF)  of  the  TTF  (TBF)  using 
Limited,  Censored  Data 

>  Censored  MLE  Approach 
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Limited  Data 


Group  LI 

Original  data  Updated  data 
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Censored  MLE  Approach 


■  Using  available  limited  data  (TBFs  and  censoring  mileage), 
“estimate”  PDF  of  TBF  using  a  censored  MLE  approach. 

■  Tail  sample  the  PDF  of  previous  step  to  “enhance”  the 
original  limited  data. 

■  Using  “enhanced”  data  from  previous  step,  “better 
estimate”  the  PDF  of  TBF  using  an  uncensored  MLE 
approach. 

■  Using  the  PDF  of  previous  step,  a  Bootstrap  approach 
estimates  statistics  of  TBF  (e.g.  distribution  of  MTBF, 
distribution  of  TBF  standard  deviation,  etc.) 
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Notation 
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Observation  /  Assumption 
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mi=Xi~fi{A,B,p,q\ 

[A  <  X.  <  B ,  and  p  >  0,  g  >  0) 

f(x,A,B,p,q)=p(p,q )  1 

[: x  —  A)p  1  (B  -  x] 
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A=  0 

B  =  45,000  miles 
p  =  3,  q  =  5 
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Observation  /  Assumption 

Beta  distribution  family  is  used  to  model  TBF. 


Beta  distributed  PDF 


A= 0,  B  =  30000 

f(x,A,B,p,q)=fi(p,qyi  (x  -  A)p~l  (B  -  x)q~l  / (B  -  A)p+q~l ,  (A  <  x  <  B,  and  p  >  0,  q  >  0) 
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MLE  Approach 

Determines  parameters  (A,  B,  p,  q )  of  “most  likely”  Beta 
distribution  using  available  data.  It  provides  Likelihood 
function  in  Bayesian  estimation. 

Censored  MLE 


#  of  recorded 
failures 


#  of  survivals 


Max  /(x,. ,  A,  B,  p,  [l  -  f(xj  ,  A,  B,  p,  <?)] 

j= i  \ 


A,B,p,q 


Beta  PDF 


Beta  CDF 


Uncensored  MLE 
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Max  FT f(xnA,B,p,q) 

A,B,p,q 
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System  (Vehicle)  Reliability 


Vehicle 

TTF 


Histogram 
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System  (Vehicle) 
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__  Monte  Carlo 
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Simulation 


For  Vehicle  : 

oo 

MTBF  =  f  R(t)dt 
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Reliability  Allocation 


Specify  system  (vehicle)  reliability 


Optimization 


V 


Determine  required  reliability  of  EACH  component 


(S 


This  optimization  problem  DOES  NOT 
have  a  unique  solution 
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Reliability  Allocation 


One  way  to  get  a  unique  solution  is  to  trade-off 
reliability  and  associated  cost 


min  Cost 


D 

—comp 


s.  t.  System  Reliability  =  R 


Target  system 
reliability 


By  varying  R'  ,  we  get  the  so  called  “Pareto  Frontier.” 
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Reliability  vs  Risk  of  Failure  (Cost) 

We  want  to  maximize  Reliability  and  simultaneously 

minimize  Risk  of  failure  (cost) 
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Reliability  -  Cost  Pareto  Front  Calculation 
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'  Definition  of  Design  Variables 

f(x,A,B,p,q)r=p{p,q)1  (x  -  A)P~\B  -  x)q~l  /{B  -  A)p+q~l ,  (A<  x  <  B,  and  p  >  0,  q  >  0) 


Beta  Distribution 


Mileage 


ju  =  MTBF 

Assume  constant  COV 
Then: 

_  U  -  A  _  (7 

U= -  <7  = - 

B-A  B-A 
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Reliability-Cost  Relation 

cost  =  COSt0  £^MTBF/MTBFo_1)  ;  por  eaci!  component 

Nc 

Cost  =  ^  [cost0 

„■  i 

e,(MTBF/MTBF0-l)(1  +  counts)]. 

tc-i 

For  system  with  Nc  components 

Unit 

COStj 

cost0 

Cost 

- ' - ' - ►  MTBF 

MTBF0  MTBFj 
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Input  Information 


Component 

Number 

Comp  No. 

Baseline  MTBF  in  hours 
(MTBFg) 

Coefficient  of 
Variation 

^ factor 

Baseline  cost 
(COStg) 

k 
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Histogram  of  System  failures 
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Reliability  Comparison  between  Repairable 
And  Non-Repair  able  System 

Reliability 
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8000 


Reliability 


System  Reliability-Cost  Pareto  Front 
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